Spin-Gap Physics, Ground State Degeneracy, and Bound States 
on the Depleted Kagome Lattice 
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We analyse the antiferromagnetic spin-i Heisenberg model on a depleted kagome lattice, where 
some bonds have been reduced to exchange integral Ji <C J\. The fully depleted system consists 
of ID chains, each with a doubly degenerate singlet-pair ground state and a spectral gap (like the 
Majumdar-Ghosh model). There are localised and itinerant low-energy excitations. The modes from 
the lowest branch of excitations are incapable of lifting the 2D system's ground state degeneracy at 
finite J2/J1 <C 1. Low-energy excitations of the 2D system are dominated by coherently propagating 
bound states of the ID excitations. 
PACS No: 75.10.Jm 



INTRODUCTION. 

In recent years, interest in systems exhibiting a spin 
gap in the excitation spectrum has been fuelled by the 
appearance of spin-gap phenomena in various situations, 
including frustrated spin systems jl). The analysis of 
toy models of spin-gap physics is noticeably easier in one 
spatial dimension, so considerable attention has been de- 
voted to this line of work. 

Despite the well known distinction between integer and 
half-integer spin excitation spectra, some spin-^ systems 
with intermediate interaction ranges have a spin gap: the 
Majumdar-Ghosh (MG) model [^|^] and several related 
models ||,[| exhibit a rich variety of behaviour. The 
study of the low-lying excitations of the MG model has 
proved a non-trivial task J6[-^0|, and we need a spin- 
gapped model analogous to the MG one, but in which 
the low-lying excitations are easily tractable. We present 
such a model, obtained by removing certain bonds from 
the spin-i Heisenberg model on a kagome lattice (cf. the 
approach of [jy] to the full 2D kagome antiferromagnet). 

Under this depletion, the lattice decouples into a set 
of one-dimensional chains, which have exact singlet-pair 
ground states (like the MG model) and whose low-lying 
excitation spectrum can be obtained analytically. The 
ground state of each chain is doubly degenerate, and con- 
sequently the ground state of the whole system has an 
degeneracy exponential in the number of chains. 

We then reinstate the removed bonds at a strength 
Ji "C Ji, where Ji is the on-chain exchange integral. 
This introduces a branch of domain wall-domain wall 
bound states which propagate coherently on the 2D lat- 
tice, and have lower energy than the on-chain modes. 
Furthermore, we show that the ground state degeneracy 
obtained for J2/J1 = is not lifted by the modes in the 
lowest-lying band for finite J2/J1 -C 1. 

THE MODEL. 

The Hamiltonian of our model is given by 



where the Si are vectors of spin-i operators. It represents 
the system shown in Fig. pi where the set A contains the 
bonds on the edges of the grey triangles. These are the 
bonds that have been weakened to strength J2/J1 with 
respect to the remainder. 




FIG. 1. The depleted kagome lattice. The grey triangles 
of bonds have been weakened to exchange integral J2; the 
remainder have exchange integral Ji 3> J2. 

THE ONE-DIMENSIONAL CHAINS. 

Ground States. In the limit J2/J1 = 0, the 2D 
system decouples into a set of ID chains. Each chain 
has doubly degenerate ground states made up of nearest- 
neighbour singlet pairs. One ground state (called the 
'red' ground state for convenience) is shown in Fig. |^. 
The other ('blue') ground state is obtained by a mirror 
reflection of the red state about the grey bar in the figure. 
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FIG. 2. The red ground state of the chain. Reflection about 
the grey bar yields the degenerate blue ground state. 

These are clearly eigenstates of the on-chain part of 
the Hamiltonian (|l|), since each bond is either occupied 
by a singlet, or part of a triangle with a singlet on one 
side; see the work of Klein |12| and Affleck et al. |13| on 
the relationship between total-spin projection operators 
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and valence bond crystal states. In the former case, an 
elementary calculation gives 
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whereas the contributions of the other two sides of the 
triangle cancel out: 



(Si • S3 + S2 • S3) 



1 



_(itix)-urx)) 



= 0, 



where X— f, I- Notice that the ground states of Fig. g 
are orthogonal to each other only in the limit when the 
number of spins on the chain tends to infinity, since 
( blue I red) ~ 2~ N , where 3N is the number of spins. 

The bonds at the top of the chain are occupied by 
singlet pairs in both ground states. If we ignore them 
for the time being (see below), we are left with a model 
rather like a Majumdar- Ghosh chain, except that only 
alternate sites have next-nearest-neighbour interactions. 

This model, obtained by removing the top bonds from 
our depleted kagome chain, is called the "delta chain" , 
or "sawtooth lattice" , and has been analysed by several 
authors. Early work Jl^l6| was numerical, and estab- 
lished the existence of an energy gap; later studies [ p7|Jl8| 
concentrated on analytic work (mainly variational calcu- 
lations), modelling the low- lying excitations as domain 
walls of a few spins. We summarise the domain- wall 
physics of a single chain here, in a notation suitable to 
the remainder of the letter. 

Low-Energy Excitations. It is easy to show that 
a localised excitation introduced into a red ground state 
will not propagate coherently. To do this, we consider a 
group of four spins on the chain (1A, IB, 2A, and 2B — 
see Fig. ||), and rewrite the local Hamiltonian in a basis 
consisting of two-spin eigenstates. Namely, 
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where the subscript n refers to the bond between the 
spins described. In this case, we consider bonds 1 and 2. 

Initially, the state of these bonds will be given by 
\S) 1 \S) 2 . Henceforth we shall drop the subscripts 1 and 
2, leaving them implicit in the order of the kets. Now let 
us consider introducing a local excitation by modifying 
the state of bond 1. It is easy to show that 
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while \S) \S) is an eigenstate of H. Since there is no 
term on the right hand side of the form \ S) \X), coherent 



propagation of such a localised excitation cannot occur. 
These considerations lead us to conclude that the low- 
energy excitations of this chain are of the domain-wall 
type. If we begin with an initially red chain, we may 
introduce a region of blue ground state into it, bounded 
by two domain walls of distinct natures (see Fig. ||). 

The stationary domain wall (boxed in the figure) is 
flanked by singlets, and hence remains spatially localised. 
The mobile domain wall (circled) feels the effect of the 
on-chain Hamiltonian, and so acquires a dispersion. 
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FIG. 3. The two types of domain wall on the chain: sta- 
tionary (boxed) and mobile (circled). 

We take our basis states to be \j), representing a mo- 
bile domain wall at site j on the chain. The calculation 
of the dispersion relation is modified by the fact that 
these basis states are not mutually orthogonal. Conse- 
quently, an expansion of our eigenstate in these basis 
states, \ip) = J2j c j \ yields an equation for the vector 
Cj which reads He = EAc, where fiy = (i|-ff c hains|i) an d 
■A43 = The energy spectrum is therefore obtained 

from the eigenvalues of the effective Hamiltonian A X H. 

Hence we find that the dispersion relation for the mo- 
bile wall is given by 
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where the ground state is taken to have zero energy. This 
result is independent of the correlation (singlet or triplet) 
between the mobile and stationary domain walls, so this 
branch contains magnetic and non-magnetic modes. 

Local Modes. As well as these domain wall modes, 
the spectrum contains modes in which the spin pairs on 
the top horizontal bonds are excited into a triplet state, 
whose role must be considered. However, compared to 
the gap of the domain wall band (w Ji/4), the singlet- 
triplet gap on these bonds (= Ji) is rather large, so we 
propose to ignore excitations of the top bond singlets 
when considering low-energy processes on the chain. 

EXCITATIONS OF THE 2D MODEL. 

Pair-Hopping Amplitude. Analysing the effect of 
a =/= on the excitations of the one-dimensional chains 
leads naturally to the consideration of tunnelling pro- 
cesses; however, it is clear that a single mobile wall can- 
not tunnel by itself, as this would violate the restriction 
that the chain states be unchanged at large distance. 

However, there exists an amplitude for pair hopping of 
a stationary and a mobile domain wall. To describe it, 
let us introduce a set of basis states more general than 
those used in the ID analysis. In our new notation, \i, j). 
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represents a state of chain 1 with a stationary domain 
wall at site i, and a mobile domain wall at site j. Our 
previous basis states \j) are now represented by \0,j) 1 . 

Let us deal with the case where both chains are in the 
red ground state at large distance. There is a term of 
perturbation theory, first order in J2/J1, that connects 
the states \i,j} 1 |red) 2 and (red^ | A; J 7'} 2 . This term is 
represented diagrammatically in Fig. 
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FIG. 4. A pictorial representation of the steps in the in- 
ter-chain tunnelling process. Those triangles through which 
the grey bar passes have bonds with exchange integral J2. 

Using the standard formulae of Rayleigh-Schrodinger 
perturbation theory, and approximating the energy of the 
local mode by Ji , we may obtain the amplitude 



Aionrefl — ( ^) 



e((i-j)(k-j)), (5) 



where Q(x) is the unit step function. The precise value 
of the prefactor depends on the value we use for the local 
mode gap, but the linearity in J2 is a basic feature. 

This amplitude corresponds to the hopping of a mo- 
bile domain wall from site j on chain 1 to the same site 
on chain 2, while the stationary wall hops from site i 
on chain 1 to a (potentially different) site k on chain 
2. The step function guarantees that the mobile wall is 
on the same side of the stationary wall after the hop- 
ping as it was before. This treatment applies equally to 
the case where the two ground states are blue. Both of 
these cases will be called 'non-reflective hopping', since 
the sign of k is preserved in the hopping process. The 
alternative situation, 'reflective hopping', will be taken 
up in a forthcoming publication |l9| ; we do not expect 
our qualitative results to be altered by the introduction 
of reflective hopping (see concluding discussion). 

Bound States (Non- Reflective Hopping). To 
analyse the non-reflective case, we formulate an effective 
model by introducing creation and annihilation opera- 
tors for the mobile (cVjC^) and stationary (d ni ,d ny i) 
domain walls; here, n and i are chain and site indices re- 
spectively. It is clear that the domain walls must obey an 



exclusion principle, but this still leaves a choice of statis- 
tics: fermions or hard-core bosons? In the dilute limit, 
the exchange statistics are not relevant to the calculation; 
conventionally, we choose the operators {c,c* ,d,d'} to 
obey Fermi statistics. 

In these terms, the effective Hamiltonian is 



ffeff =t 
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(cVc„, i+ i + h.c. 
+ (e - M) c U Cn >* + ^ Y d l,i d n,i 

n,i n,i 
+ J2 {Kk hn 4 t ,i d i,k d n-l,jCn-l,i+ ]x - c ) ( 6 ) 

where t = X J\, z = lf«/i, and Vijk is just the amplitude 
^4-nonrefl from Eq. (||), with i being the position of the 
mobile wall, j that of the stationary wall before hopping, 
and k that of the stationary wall after hopping. 

Because of the translational invariance in the x (paral- 
lel to the chains) and y (transverse to the chains) direc- 
tions, we may Fourier transform and work in momentum 
space. The action is given by 

S = Y c lph - E - It COSp + H) C^ph 

+ £ r(0) SK» U , C t 3 rh d L S h d ^lhC^ph, (7) 
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where UJ4 — uj\ + U2 — s — p + q 



r, and ct^, c up h, 

wp /ii d^ph now represent Grassmann fields. The Wi are 
fermionic Matsubara frequencies, p, g, r, s are momenta 
along the chains, and h is the transverse momentum. 
Recall that, in the case we consider, all momenta are 
restricted to positive values. The bare vertex part, T^, 
is given by 



p(o) 
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(f(q)f(s) + f(-q)f(s)), 



(8) 



where the function f(q) is the Fourier transform of the 
hopping amplitude's spatial dependence, given in Eq. (|^). 

The formation of bound states is signalled by the ap- 
pearance of divergences in the dressed vertex part, T. 
Being unable to solve the full set of Schwinger-Dyson 
equations, we make a simple Ansatz: we assume the self- 
energies of the two propagators to be zero. Our defence 
is that the gap in the unperturbed model is quite large. 

The remaining bound state condition can be solved 
analytically only near two lines in the Brillouin zone, h — 
±77- Setting h = \ — 5, we find that 



J 2 S I 2 cos 2 P 
m \ 



4 cos P + 3 
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where P is the on-chain momentum, and ^e 2 is the dis- 
tance that the bound state lies below the band of un- 
perturbed modes, i.e. Abound = 52^1 ~~ \^ ■ Similar 
behaviour is found near the line h = — \. We tentatively 
conclude that the branch of bound states exists through- 
out the Brillouin zone, touching the bottom of the con- 
tinuum band tangentially along the lines h = ±5. 

GROUND STATES OF THE 2D MODEL. 

In the limit Jij J\ = 0, the double degeneracy of each 
chain in the system leads to a total degeneracy g = 2 L , 
where L is the number of chains. 

It is noteworthy that we have been unable to find a pro- 
cess (within the manifold of domain wall states) that lifts 
this ground state degeneracy in the a^O case. Any such 
process consists of a finite number of intcr-chain hopping 
events. Each such hop between two chains of the same 
ground state has a counterpart hop between two chains 
of different ground states, the only difference being that 
in the latter case the left-to-right order of the stationary 
and mobile domain walls is reversed. Due to a combina- 
tion of the spatially extended nature of the domain wall 
states (see above) and the imposition of periodic bound- 
ary conditions, such a reversal cannot make a difference 
to the amplitude of the process. Hence all perturbation 
theory corrections to the energy of a given ground state 
of the two-dimensional lattice have counterparts of equal 
magnitude and sign for any other ground state. 

Although no explicit calculation has been performed, 
the degeneracy is almost certainly lifted by perturbative 
processes involving the excitation of local modes on the 
'top' horizontal bonds of the chains. Nonetheless, it is 
interesting to note that, because of the effective screening 
of the chains from one another due to the singlets on 
the top bonds, the modes in our low-energy manifold are 
incapable of lifting the ground state degeneracy of the 
a = case. 

CONCLUSION. 

We have presented a quasi-lD chain with nearest- 
neighbour antiferromagnetic exchange and a spin-gapped 
excitation spectrum that can be treated analytically in 
the low-energy limit. This chain was obtained by a nat- 
ural depletion of the 2D kagome antiferromagnet, and as 
well as itinerant and localised domain walls (obtained in 
|p7| , ^8| by similar methods to ours) it has local modes 
with a gap larger than that of the itinerant excitations. 

Recoupling these chains weakly, we have obtained a 
pair-hopping amplitude of the itinerant and localised ID 
modes that arises naturally from the nature of the ID 
ground states. Furthermore, we have shown that the 
ground state degeneracy of the 2D system (which is expo- 
nential in the number of chains) is not lifted by processes 
in the lowest branch for finite J2/J1 <C 1; this effect is 
due to symmetry properties of the ground states and the 
nature of the interchain coupling. As regards the exci- 
tations, the pair hopping amplitude leads to a band of 
bound states of the ID domain walls: these bound states 



propagate coherently on the 2D lattice, and have lower 
energy than the on-chain excitations. 

While our model is not proposed as a starting point 
for treating the isotropic 2D kagome antiferromagnet, it 
has several novel and interesting features, notably the 
inability of the domain wall modes to lift the ground state 
degeneracy p0| and the formation of a band of coherently 
propagating 2D states for arbitrarily small J2/J1 7^ 0. 

Finally, one might ask whether the consideration of 
reflective rather than non-reflective hopping will change 
our results. Qualitatively, one may argue that it will 
not: since bound states already form in the non-reflective 
system, the existence of reflective hopping will act only to 
strengthen the tendency to bind, perhaps modifying the 
prefactor in (||) . These topics, along with a semiclassical 
method of treating the non-reflective case, will be taken 
up in a forthcoming publication [fl9|| , 
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